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Abstract 

■ 

ON ', Motivated by the recent observation of double charm quark pair production by 

the Belle Collaboration, we calculate the complete 0(a 2 s ) inclusive production cross 
■ sections for r] c , J/ip, and Xcj(J=0, 1, 2) plus cc in e + e~ annihilation through a vir- 

tual photon. We consider both color-singlet and color-octet contributions, and give 
the analytical expressions for these cross sections. The complete color-singlet calcula- 
£SJ ■ tions are compared with the approximate fragmentation calculations as functions of the 

center-of-mass energy yfs. We find that most of the fragmentation results substantially 
overestimate the cross sections (e.g. by a factor of ~4 for Xci and Xc2) at the Belle and 
BaBar energy ^fs = 10.6GeV. The fragmentation results become a good approxima- 
. tion only when y/s is higher than about lOOGeV. We further calculate the color-octet 

Qh! contributions to these cross sections with analytical expressions. We find that while 

O 1 the color-octet contribution to J/ ip inclusive production via double charm is negligible 

(only about 3%), the color-octet contributions to Xci and Xc2 can be significant. 
PACS number(s): 12.40.Nn, 13.85.Ni, 14.40.Gx 
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: 1 Introduction 



Charnonium is one of the simplest quark-antiquark composite particles. Charmonium physics 
has played an important role in the study of Quantum Chromodynamics (QCD) both pertur- 
batively and nonperturbatively, since the first charmonium state J/ip was discovered in 1974. 
During the past decade, the study of charmonium has become more interesting because of 
the large difference between the predictions of the color-singlet model and the observations 
of J/i/j and production at several experimental facilities e.g. at the Fermilab Tevatron p^. 

The newly developed nonrelativistic QCD (NRQCD) factorization formalism [2j allows 
the infrared safe calculation of inclusive heavy quarkonium production and decay rates. In 
the NRQCD production mechanism, a heavy quark-antiquark pair can be produced at short 
distances in a conventional color-singlet or a color-octet state, and then evolves into an 
observed quarkonium nonperturbatively. With this color-octet mechanism, one may explain 
the Tevatron data on the surplus production of J/ip and ip' at large px, though puzzles about 
their polarizations still remain (for a recent review see jH] and references therein). 
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To further test the color octet mechanism, it is interesting to study the charmonium 
production in e + e~ annihilation. The J /if) inclusive production in e + e~ annihilation has 
been investigated within the color-singlet model [HOEE] and the color-octet model 0IH1III]. 
The angular distribution and energy distribution of color-singlet J /if) production at s/s = 
10.6 GeV have been discussed in jH]. In it is found that a clean signature of the color- 
octet mechanism may be observed in the angular distribution of J /if) production near the 
end point region. In jSj contributions of various color-octet as well as color-singlet channels 
to the J /if) production cross sections are calculated in a wide range of e + e~ collider energies. 
Moreover, the J /if) polarizations are predicted in Recently, BaBar [TTJ and Belle [12J have 
measured the direct J /if) production in continuum e + e~ annihilations at y/s = 10.6 GeV. 
The total cross section and the angular distribution seem to favor the NRQCD calculation 
over the color-singlet model ^JJ, but some issues (e.g. about the momentum distribution 
and polarization of J /if)) still remain. 

The situation has become even more complicated due to the very recent observation 
for the double cc production associated with J /if) by Belle [EI]. The measured exclusive 
cross section for e + + e~ — > J/if> + r\ c process is an order of magnitude larger than the 
theoretical value [J3], and the measured inclusive cross section for e + + e _ — ► J /if) + cc 
(~0.9 pb) |12t IT3*j is more than five times larger than NRQCD predictions which are only 
about 0.1-0.2 pbjHl El M UD] taking into account the differences in the values of the input 
parameters or methods. Among other attempts to solve the J/if>cc inclusive production 
problem [T5] ITH]. e + e~ annihilation into two photons was also studied [d HE], but the two 
photon contribution turned to be negligible at \fs = 10.6GeV, though it could prevail over 
one photon contribution at higher energies (say, \fs > 20GeV)[T7J. 

The double cc production associated with J / if (both exclusively and inclusively) is very 
puzzling and needs a better understanding for both perturbative and nonperturbative QCD. 
On the other hand, experimentally, it is not clear whether the copious (even dominant) 
double cc production will also happen for charmonium states other than J /if), such as r] c 
and Xcj(J=0, 1, 2). Among them the x& and Xc2 are more interesting since they have 
large branching fractions decaying into J /if) + 7 and might be easier to be detected. In 
fact, in Ref. [12] the inclusive production for Xci and Xc2 was searched for with the available 
integrated luminosity of about 30 fb _1 at Belle. As more data are collected in the near 
future at B factories we hope that more accurate measurements for the P-wave and other 
S-wave charmonium states will be possible. These measurements will be helpful to clarify 
the problems associated with J /if) double cc production. 

On the theoretical side, the calculations for inclusive S-wave and P-wave charmonia 
production via double cc are necessary in the framework of NRQCD, including both the 
color-singlet and color-octet contributions. When we know the differences between NRQCD 
predictions and experimental data, we will have to further consider other mechanisms and 
methods in QCD to explain the differences. 

This paper is organized as follows. In Section 2, we will calculate the complete 0{o? s ) 
color-singlet inclusive production cross sections for r\ c and Xcj{ J=0, 1,2) (along with J /if)) 
via double cc in e + e~ annihilation through a virtual photon. Then we will compare the 
complete calculation with the calculation obtained in the charm quark fragmentation limit, 
and give their ratio as functions of the center-of-mass energies and determine the energy 
scales at which fragmentation approximations become reliable. In section 3, we will further 
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estimate the color-octet contributions to J /if) and Xcj inclusive cross sections via double 
charm. Finally, we summarize our results in section 4. 



2 Color-singlet contribution to charmonium produc- 
tion via double cc in e + e~ annihilation 

The quarkonium can be described in term of Fock states superposition within the NRQCD 
framework as follows 

| ^ Q > = 0(1) | QQCS?} > +0(v)QQ( 3 Pf ) }g > 

+0(v 2 )QQ[ 1 S^}g > +0(v 2 )QQ[ 3 S[ 1 %g >+..., (1) 

I Xcj > = 0(1) | QQ[ 3 PP] > +0(v)QQ[ 3 s[% > 

+0(v 2 )QQfpf%g >+..., (2) 

where the superscript 1 or 8 labels the color configuration of the Fock Components. 

Following the nonrelativistic QCD (NRQCD) factorization formalism, the scattering am- 
plitude of the process e~(pi) + e + (p 2 ) —>■ 7* — > cc( 2S ' +1 Lj' 8a ' > )(p) + c(p c ) + c(p e ) in Fig. [T] and 
Fig. 121 is given by 

A(e~( Pl ) + e + {p 2 ) - ccC S+1 L ( j M) )(p) + c(p c ) + c(p £ )) = ^EEE 

L Z S Z s\si jk 

x ( Sl ;s 2 I SS Z )(LL Z ;SS Z | JJ z )(3j;3k | 1,8a) 

' A{e-{pi) + e+(p 2 ) -> Cj (§; si) + c fc (§; s 2 ) + Q(f ; s 3 ) + Q(f ; s 4 )) (L = 5), 

e* a (L z )A a (e-( Pl ) + e + (p 2 ) -> c,(§; Sl ) + c fc (|; s 2 ) + c,(f ; s 3 ) + c,(f ; s 4 )) (L = P). 

(3) 

where cc( 25+1 Lj 1 ' 8a ' 1 ) is the intermediate cc pair produced at short distance, which subse- 
quently evolves into a specific charmonium state at long distance, A a is the derivative of 
the amplitude with respect to the relative momentum between the quark and anti-quark in 
the bound state. For the case of color-singlet state, the coefficient Cl can be related to the 
origin of the radial wave function (or its derivative) of the bound state as 

C s = ^-\ R s (0) | 2 , C P = ^-\ R' P (0) | 2 . (4) 



x 



47r ~ Arc 

The spin projection operator can be defined as[P3] 

>V , \-,V 

51 S 2 



p ss z {p;q) = J2(si; s 2 \SS z )v(- + q; si)u(- - q; s 2 ). (5) 



We list the spin projection operators and their derivatives with respect to the relative 
momentum, which we will use in the calculations, as 

P 00 (p,0) = ^= 7 5(i> + 2m c ), (6) 
3 



Pis,(p,0) = -jj=4(S M )(j> + 2m c ), (7) 

A\(P,0) = i - ? i—[ 7 a 4*(Sz)(fi + 2rn c )-(fi-2m c ) /(^) 7 1- (8) 
For P-wave states we need further relations to reduce the polarizations 

E e™(L z )e*P(S z )(lL z ; 1S Z | J = J z = 0) = ^(-^ + *§£) 

E ^(L z )e^(S z )(lL z ; l£ z | J = U z ) = -^-e a ^p K e\{J z ) (9) 
L Z S Z V2M 

E e* Q (L z )e^(^)(lL z ; 1S Z \ J = 1J Z ) = e*^(J z ) 

LzSz 

where M is the mass of the charmonium, which equals to 2m c in the nonrelativistic approx- 
imation. 

The calculation of cross sections for e~ + e + — > 7* — > charmonium +cc is straightforward. 
Using the definition in Ref. [2J we get the differential cross section as follows 

da(e + + e~ — ► 7* — > charmonium + cc) 4C < La 2 a: 2 . Q! (- 2 )n 

di = 81m c + — )■ ^ 

where L=S for S-wave, L=P for P-wave and z = 2E J / 1 p/^fs. The functions S(z) and a(z) 
for different charmonium states are given in the Appendix. 

With Eq. (fTUjl we can evaluate the inclusive cross sections for r) c , J/ip and Xcj- The input 
parameters used in the numerical calculations are[2*U] 

m e = 0, m c = 1.5GeV, a s {2m c ) = 0.26, a = 1/137, (11) 

I R s (0) | 2 = 0.81 GeV 3 , I Rp(0)' | 2 = 0.075 GeV 5 . (12) 
Now we give the numerical results at the Belle and BaBar energy ^fs = 10.6 GeV. 

(x(e + + e~ -> 7* -> r? c + cc) = 58.7 fb (13) 

a(e + + e~ -> 7 * -> J/V> + cc) = 148 fb (14) 

a(e + + e~ -> 7 * -> Xco + cc) = 48.8 fb (15) 

a(e + + e" -> 7* -> + cc) = 13.5 fb (16) 

a(e + + e~ -> 7 * -»■ x c2 + cc) = 6.30 fb (17) 

The J/i/j production rate is in agreement with other references [HI El HI] after taking into 
account the differences in the values of the input parameters. In the z ^> 5 limit, where 5 is 
defined as 4m c / yfs, the approximate fragmentation results will be equivalent to the complete 
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calculations. This is another check for the validity of the complete calculation. Here the 
fragmentation cross sections are written as 



<Tf rag (e + e — > 7* — > charmonium + cc) = 

2a(e + + e~ -> cc) / T> c ^ charmonium (z)dz, (18) 

where P(-z) are the charm quark fragmentation functions into S-wave [2~T] or P-wave[22] 
charmonia. 

The cross sections obtained in the complete calculation and in the fragmentation approx- 
imation as functions of the center-of-mass energies are plotted in Fig. EE All these cross 
sections are in units of a cc = a(e + + e~ — > 7* — > cc), the cross section for e + e~ annihilat- 
ing into the cc quark pair, times 10~ 4 . One can find that the cross sections in complete 
calculations and fragmentation approximations (all in units of the cross section for e + e~ 
annihilating to the cc pair) are proportional to the fragmentation probabilities for the charm 
quark fragmentating into charmonia when the (5 C 1 limit is valid. This is just what the 
fragmentation approach describes. The results in these figures show that except for x c0 , the 
differences between fragmentation results and complete calculations are large at low ener- 
gies. At the Belle and BaBar energy ^/s=10.6GeV, the ratios of complete calculations to 
fragmentation calculations are 

+ 6 ~ ^ 7 * charmomum + cc ") = o.28, 0.58, 0.25, 0.25 (19) 



<jf rag (e + + e — > 7* — ► charmonium + cc) 

for n c , J/ip, Xci, and Xc2 respectively. As the center-of-mass energy increases, the ratios 
of complete calculations to fragmentation results increase and can reach to 90% when the 
center-of-mass energy is over 100 GeV. Moreover, the cross sections are rather sensitive to the 
input parameters. If we choose a = 1/134, a s = 0.28, and m c = 1.48GeV at y/s = 10.6GeV, 
the cross sections for r/ c , J/ip, and \cj (J=0,l,2) become 77.0fb, 192fb, 64.2fb, 18.3fb, 8.48fb 
respectively. 

3 Color-octet contribution to J/ip and \cJ production 
via double cc in e + e~ annihilation 

We next estimate the color-octet contribution to J/ ip and Xcj production via double cc in 
e + e~ annihilation. The Feynman diagrams are showed in Fig ^ and Fig|21 

In Fig. ^ the charmonium comes from the color-octet mediate states cc( 2S+1 L^) by 
emitting soft gluons at long distances. Here the color-octet contribution can be obtained 
from the corresponding color-singlet contribution divided by a factor of ^§k^w+TZJ^ ■ The 
matrix elements < ( 2S+1 Lj) > can be extracted from the Tevatron data for J/ifi and 
Xcj production (see Ref. [21] for detailed discussions). Accordingly, with some unavoidable 
uncertainties we set them as 

< Of ( 3 Si) >= 1.16 GeV 3 , (20) 
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< 0% cl ( 3 P 1 ) >= 0.32 GeV 5 , 



(21) 



< C|( 3 5i) >= 1.06 x 10~ 2 GeV 3 , (22) 

< Of( 3 P ) > /m 2 c = 1.0 x 10~ 2 GeV 3 , (23) 

< 0%( 3 Pj) >= (2 J + 1) < Of ( 3 P ) >, (24) 

< Ogf cl ( 3 5i) >= 1.0 x 10" 2 GeV 3 . (25) 

With these values of the matrix elements, the color-octet contributions to J/ip and Xcj 
in Fig. ^ are about at least two orders of magnitude smaller than the color-singlet contri- 
butions, and therefore are negligible. In Fig. E] the color-octet contributions come from four 
different (the upper two and the lower two) diagrams. With their contributions (including 
the interference terms), the differential cross section reads 

dogctet = I6a 2 a 2 s < Og( 3 Si) > , |2 
dz 27 m c 

where I M I 2 takes the form 



M \\ (26) 



- 2 7T (l-z)(z 2 -5 2 ) 

1 ~ 125 2 s 2 z{z-2) 2{ Z \ A + 5 2 -Az 

[3d 4 - 125 2 (z - 2) + 16(10 + z(z - 10))] + 
(z - 2) 2 [35 4 - 85 2 (3z - 4) + 32(2 + z(z - 2))] 



zV4 + 5 2 - Az + 2,/(l - z)(z 2 - 5 2 ) 

ln[ V = ]}■ (27) 

zV4 + 5 2 - Az - 2J{l-z){z 2 -5 2 ) 



The numerical results can be obtained by using the parameters given above, and are 

<T oc tet{e + e- -> J/^cc) = 4.5 fb. (28) 

cr oc tet(e + e~ Xacc) = 4.3 fb. (29) 

The color-octet contribution to J/ip is only 3% of the color-singlet cross section. For 
Xdi the color-octet contribution is significant, which is about 32% of the color-singlet cross 
section. With the approximation of heavy quark spin symmetry, the contributions of color 
octet 3 Si to Xcj (from color-octet 3 Si mediate state to color-singlet 3 Pj final state by El 
transition) satisfy the ratio 1:3:5 for J = 0, 1, 2 respectively. Their values are given by 

o"orfe*(e + e" -> Xcocc) = 1.4 fb. (30) 
&octet(e + e~ -> Xc2cc) = 7.2 fb. (31) 
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We show the angular distribution and energy distribution for Xci in Fig |S] and Fig El One 
can see that the color-octet contribution enhances the differential cross section significantly 
in the low energy (small z) region. 



4 Conclusion 

In summary, we have calculated the complete 0(a 2 ) inclusive production cross sections for r) c , 
J/ip, and Xcj(J=0, 1, 2) plus cc in e + e~ annihilation through a virtual photon. We consider 
both color-singlet and color-octet contributions, and give the analytical expressions for these 
cross sections. The complete color-singlet calculations are compared with the approximate 
fragmentation calculations as functions of the center-of-mass energy yfs. We find that most 
of the fragmentation results substantially overestimate the cross sections (e.g. by a factor 
of ~4 for Xci and \c2) at the Belle and BaBar energy ^fs = 10.6GeV. The fragmentation 
results become a good approximation only when -y/i is higher than about lOOGeV. We further 
calculated the color-octet contributions to these cross sections with analytical expressions. 
We find that while the color-octet contribution to J/ip inclusive production via double charm 
is negligible (only about 3%), the color-octet contributions to Xci and Xc2 can be significant. 
These results may serve as NRQCD predictions to compare with the experimental data 
observed or to be observed at Belle and BaBar. 
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Appendix 

In this Appendix, we give the functions of S and a which are defined in Eq. (|10p. 

Vc 3s 2 5 2 z%z-2)%z 2 -5 2 yl A + 5 2 -Az I [2 + d ){4 + d , 
+96<5 6 (64 + 225 2 + 5 4 )z - 165 2 (1920 - 8QA5 2 + 532<5 4 + 1255 6 - 25 8 ) 



z 2 



)z 3 



,-85 2 (9984 - 53125 2 + 488<5 4 + 96<5 6 - 5\ .. 
+2(6144 - 47872<5 2 + 208005 4 - 3925 6 - 1105 8 + 35 10 )z 4 
-4(6144 - 213765 2 + 42565 4 + 1125 6 + 95 8 )z 5 
+ (14336 - 513285 2 + 54725 4 + 4205 6 - 35 8 )z 6 
-4(1536 - 31685 2 + 3525 4 + 5 6 )z 7 
+8(864 - 3Q5 2 + 135 4 )z 8 - 32(112 + 115 2 )z 9 + 768^ 10 ] 
-35 2 (z - 2) 4 [8<5 6 (2 + 5 2 ) - 9Q5 6 z - 25 2 (192 - 485 2 + 85 4 - 5 6 )z 2 
+165 2 (8 + 65 2 - 5 4 )z 3 + <5 2 (192 + 40<5 2 - <5 4 )z 4 + 8(32 - 4<5 2 + 5 A ) 
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, „ zVA + S 2 - Az + 2 J {I - z){z 2 - 5 2 ) 

-8(48 + 5 2 )z 6 } In V = 

z^A + 5 2 -Az- 2^(1 - z){z 2 - 5 2 ) 



Tlf ^ [96^(4 + 0(6 + 



3s 2 (5 2 2 3 (2-2) 6 (^ 2 -<5 2 ) L V 4 + <5 2 -4z 
-965 6 (64 + 185 2 + <5 4 )^ + 16<5 2 (2688 + 6085 2 + 4285 4 + 435 6 - 25 s )z 2 
-85 2 (17664 + 32645 2 + 1845 4 - 9Q5 6 - 5 8 )z 3 
+2(6144 + 893445 2 + 77445 4 - 20245 6 - 1745 8 - 35 w )z 4 
-4(6144 + 226565 2 - 13765 4 - 5125 6 - 355 8 )z 5 
+ (14336 + 55045 2 - 51525 4 - 7325 6 - 35 8 )z 6 
-4(1536 - 17605 2 - 4165 4 + 5 6 )z 7 
+8(864 - 1965 2 + 135 4 )/ - 32(112 + 115 2 )^ 9 + 768^ 10 ] 
+35 2 (z - 2) 4 [85 6 (6 + 5 2 ) - 325 6 z - 25 2 (QA + 485 2 + 16<5 4 - 5 6 )z 2 
+ 165 2 (12 - 5 2 )(2 + 5 2 )z 3 - (1024 + 3205 2 - 885 4 - 5 A )z A 
+8(96 - 285 2 - 5 4 )z 5 + 8(16 + 5 2 )z 6 ] 

z^A + 5 2 -Az + 2^(1 - z){z 2 - 5 2 ) 
X ln zy/4 + 5 2 -Az - 2^(1 -z)(z 2 -5 2 )^' 



An _ { ^/(l-,)(, 2 -5 2 ) 



r8 U 2 



s 2 5 2 z 3 (z - 2) 6 (^ 2 - 5 2 ) 1 V 4 + 5 2 - Az 
x [-325 4 (4 + <5 2 )(48 + 225 2 + 35 4 ) 
+32<5 4 (768 + 4005 2 + 665 4 + 35 6 )z 
-165 2 (384 + 19205 2 + 5565 4 + 295 6 - 25 
+85 2 (1792 + 1285 2 - 568<5 4 - 805 6 - 5 8 )z 3 
+2(2048 - 110085 2 + 107525 4 + 31765 6 + 985 8 + 35 w )z 4 
-4(4096 - 78085 2 + 3424<5 4 + 6005 6 + 175 8 )^ 5 
+ (38912 - 206085 2 + 4544<5 4 + 5085 6 - 35 8 )z 6 
-4(13312 - 8005 2 + 1205 4 - 35 6 )z 7 + 8(4512 - 205 2 - 155 4 )z 
-32(336 - 5 2 )z 9 + 12802 10 ] 

-5 2 (z - 2) 4 [8<5 4 (48 + 225 2 + 35 4 ) - 32<5 4 (24 + 55 2 )z 
-2<5 2 (448 + 165 2 + 85 4 - 35 6 )^ 2 + 165 2 (56 - 105 2 - 55 4 )z 3 
+<5 2 (1152 + 2725 2 - 35 4 )^ 4 + 8(32 - 925 2 + 55 4 )^ 5 - 56(16 + 5 2 )z 
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zV A + 5 2 -Az + 2^(1 -z)(z 2 -S 2 ) 

+512z m — — — \. 

z^A + 5 2 -Az- 2y/(l - z)(z 2 - 5 2 ) 



4tt / (1 - z)(z 2 -T 2 ) 

s 2 S 2 z 3 {z-2) e {z 2 -5 2 ) { Z V 4 + 5 2 - Az 



x [325 4 (4 + 5 2 ) (16 + 25 2 + 3S A ) - 325 4 (256 + 485 2 + 225 4 + 35 6 )z 

+165 2 (1152 + 10245 2 - 1405 4 - 53<5 6 - 25V 

-85 2 (5376 + 1285 2 - 15765 4 - 2405 6 - 5 8 )z 3 

+2(2048 - 7685 2 - 199685 4 - 69685 6 - 3505 8 - 35 w )z A 

-4(4096 - 200965 2 - 111685 4 - 12085 6 - 435 8 )/ 

+ (38912 - 753925 2 - 169605 4 - 9965 6 - 35 8 )z 6 

-4(13312 - 63045 2 - 8725 4 - 35 6 )z 7 + 8(4512 - 500<5 2 - 155 4 )/ 

-32(336 - 5 2 ) + 1280z 10 ] 

+5 2 (z - 2) 4 [85 4 (16 + 25 2 + 35 4 ) - 325 4 (8 - 5 2 )z 

-25 2 (320 - 2725 2 + 64<5 4 - 35 6 )^ 2 + 165 2 (40 - 545 2 - 55 4 )z 3 

-(1024 - 7205 4 - 35 6 )^ 4 + 8(96 - 365 2 - 55 4 )z 5 



„ R 7 zy/4 + 5 2 - 4z + 2 J(l - z)(z 2 - 5 2 ) 

+8(80 + 7S 2 )z 6 - 512^ 7 ] In V = =}. 

zy/A + 5 2 - 4z - 2^(1 - z){z 2 - 5 2 ) 



8tt 



9sWz 5 (z-2) & (z 2 -5 2 ) 



{-4zy/(l - z){z 2 - 5 2 )(A + 5 2 - Az) 



x [23045 10 - 11525 8 (26 + 55 2 )z + 1925 6 (640 + 4645 2 + 355 4 )^ 2 
+965 4 (1152 - 48165 2 - 11365 4 - 435 6 )z 3 
+ 165 2 (4608 - 330245 2 + 44752<5 4 + 43605 6 + 755 8 )/ 
-8<5 2 (21504 - 1233925 2 + 78448<5 4 + 28845 6 - A55 8 )z 5 
-4(12288 - 1566725 2 + 2442245 4 - 781285 6 - 5125 s + 215 10 )/ 
-2(24576 + 5498885 2 - 3560965 4 + 417445 6 + 805 s - 95 10 )z 7 
-8(4608 - 93952<5 2 + 45728<5 4 - 12065 6 + 275 8 )z 8 
+ (487424 - 2083845 2 + 1194245 4 + 6965 6 - 95 8 )z 9 
-4(155904 + 41605 2 + 52165 4 - 215 6 )z w 

+ 16(22976 + 14805 2 + 855 4 )^ 41 - 480(232 + 115 2 )^ 12 + 15104^ 13 ] 

+35 2 {z - 2) 4 [-192<5 10 + 965 8 (26 + 35 2 )z - 645 6 (160 + 755 2 + 35 4 )z 

-165 4 (576 - 16645 2 - 1835 4 - 25 6 )z 3 

-45 2 (1536 - 46085 2 + 40165 4 + 1525 6 + 55 8 )^ 4 

+25 2 (11264 - 234245 2 - 1605 4 + 106<5 6 + 35 8 )z 5 

+4(2048 - 42245 2 + 99525 4 + 248<5 6 - 275 8 )z 6 

-(20480 - 22528<5 2 + 53125 4 - 3685 6 + 35 8 )^ 7 

+4(4096 - 64965 2 - 600<5 4 + 17<5 6 )/ - 16(320 - 4725 2 + 75 4 )/ 

„ in zy/A + 5 2 - Az + 2 J (I - z)(z 2 - 5 2 ) 

+32(48 + 5 2 )z 10 ] In V / = =}. 

zy/A + 5 2 -4z- 2^(1 - z){z 2 - 5 2 ) 



8tt 



9s 3 5 i z 5 (z-2) 8 (z 2 -5 2 ) 



{Azy/(1 - z)(z 2 - <5 2 )(4 + 5 2 - Az) 
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x [23045 10 - 57605 8 (6 + 5 2 )z + 1925 6 (896 + 4245 2 + 355 4 )/ 
+965 4 (384 - 45285 2 - 9045 4 - 435 6 )/ 
+ 165 2 (1536 + 38405 2 + 235365 4 + 29925 6 + 755 8 )/ 
-85 2 (52224 + 87045 2 + 32805 4 + 19245 6 - 455 8 )/ 
+4(12288 + 706565 2 - 512005 4 - 342245 6 - 2325 s - 215 10 )/ 
+2(24576 + 3368965 2 + 1338885 4 + 539045 6 + 3765 s + 95 10 )/ 
+ 16(2304 - 627205 2 - 112805 4 - 21915 6 - 305 s )/ 
-(487424 - 6056965 2 - 613125 4 - 70165 6 - 95 s )z 9 
+4(155904 - 407685 2 - 25605 4 - 215 6 )/ 

-16(22976 - 5045 2 + 855 4 )/ 1 + 480(232 + 115 2 )/ 2 - 15104/ 3 ] 

+35 2 (,2 - 2) 4 [1925 10 - 2885 8 (10 + 5 2 )z + 645 6 (224 + 595 2 + 35 4 )/ 

+165 4 (192 - 12485 2 - 1215 4 - 25 6 )/ 

+45 2 (512 - 122885 2 + 23845 4 - 565 6 + 55 s )/ 

+25 2 (3072 + 359685 2 - 1605 4 + 505 6 - 35 s )/ 

-4(2048 - 49925 2 + 82245 4 - 4085 6 - 235 s )/ 

+ (12288 - 512005 2 - 30085 4 - 14565 6 - 35 8 )/ 

-4(2048 - 89925 2 - 12245 4 - 175 6 )/ + 16(192 - 6165 2 - 75 4 )/ 



,, z/4 + 5 2 - Az + 2/(1 - z)(z 2 - 5 2 ) 

+32(16 + 5 2 )/ ] In V = 

r/4 + 5 2 - 4z - 2^(1 - z){z 2 - 5 2 ) 



-11525 6 (192 + 2085 2 + 625 4 + 55 6 )^ 
+1925 4 (3072 + 64005 2 + 35685 4 + 6685 6 + 355 8 )/ 
-965 4 (26624 + 278085 2 + 99925 4 + 12765 6 + 435 8 )/ 
+165 2 (36864 + 2772485 2 + 1952965 4 + 504645 6 + 44065 s + 755 10 )/ 
-85 2 (258048 + 5219845 2 + 3026245 4 + 578005 6 + 26725 s - 455 10 )/ 
-4(98304 - 7536645 2 - 5649925 4 - 3100485 6 - 377365 s - 1725 10 + 215 12 / 
+2(983040 - 6594565 2 + 844805 4 - 1030085 6 - 90005 s - 2205 10 + 95 1 
-(4784128 + 23306245 2 + 15285765 4 + 1208005 6 + 3965 s + 1175 10 )/ 
+ (6914048 + 39280645 2 + 11377925 4 + 745445 s + 19005 8 - 95 10 )/ 
-2(3100672 + 12943365 2 + 2006725 4 + 80365 6 - 95 8 )/ 
+8(443392 + 1169925 2 + 80485 4 + 355 6 )/ 1 

-64(20288 + 28085 2 + 515 4 )/ 2 + 512(544 + 335 2 )/ 3 - 28672/ 4 ] 
-35 2 (^ - 2) 4 [-1925 8 (3 + 5 2 ) + 965 6 (48 + 285 2 + 35 4 )^ 
-165 4 (768 + 8085 2 + 2175 4 + 125 6 )/ 
+165 4 (1600 + 6525 2 + 1055 4 + 25 6 )/ 



^ 12 U 7 
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+45 2 (7168 - 43525 2 - 3605 4 - 595 6 - 55 s )/ 

-25 2 (24576 - 39685 2 + 10245 4 - 645 6 - 35 s )/ 

+5 2 (17408 - 72965 2 + 1365 4 - 515 6 )/ 

-(8192 - 128005 2 - 85765 4 - 3005 6 + 35 s )/ 

+2(8192 - 66565 2 - 13285 4 + 175 6 )/ - 128(80 - 105 2 + 5 4 )/ 

9 in z/4 + 5 2 - Az + 2/(1 - z)(z 2 - 5 2 ) 

+128f24 + 5£ 2 U 10 ]ln V = ==}. (38) 

zy/A + 5 2 -Az- 2^(1 - z){z 2 - 5 2 ) 



-11525 6 (64 + 805 2 + 425 4 + 55 6 / 

+ 1925 4 (1024 + 24325 2 + 13605 4 + 4045 6 + 355 s )/ 

-965 4 (8192 + 118725 2 + 36405 4 + 6525 6 + 435 s )/ 

+165 2 (110592 + 586245 2 + 713285 4 + 128645 6 + 15225 s + 755 10 )/ 

-85 2 (724992 - 2455045 2 - 210245 4 - 725 6 + 1365 s - 455 10 )/ 

+4(98304 + 13926405 2 - 19002885 4 - 3493445 s - 206485 s - 17565 10 - 215 12 )/ 

-2(983040 - 8560645 2 - 50785285 4 - 7043525 6 - 377365 s - 7245 10 - 95 12 )/ 

+ (4784128 - 73523205 2 - 74129925 4 - 7607365 6 - 204525 s - 4475 10 )/ 

-(6914048 - 61972485 2 - 32254725 4 - 2025765 6 - 47725 s - 95 10 )/ 

+2(3100672 - 12431365 2 - 3768645 4 - 152605 6 - 95 s )/ 

-8(443392 - 473605 2 - 101285 4 + 355 6 )/ 1 + 64(20288 + 4725 2 + 515 4 )/ 2 

-512(544 + 335 2 )/ 3 + 28672/ 4 ] 

+35 2 / - 2) 4 [-1925 s (l + 5 2 ) + 965 6 (16 + 125 2 + 35 4 / 

-165 4 (256 + 3445 2 + 595 4 + 125 6 )/ + 165 4 (448 + 4045 2 - 95 4 + 25 6 )/ 

+45 2 (5120 + 17925 2 - 8565 4 + 1355 6 - 55 s )z 4 

-25 2 (16384 + 60165 2 - 10885 4 + 245 6 - 35 s )/ 

+(32768 - 256005 2 - 70405 4 - 18645 6 - 575 s )/ 

-(57344 - 814085 2 - 119045 4 - 8845 6 - 35 s )/ 

+2(16384 - 276485 2 - 23845 4 - 175 6 )/ - 128(48 - 1185 2 - 5 4 )/ 

„ in zVA + 5 2 -Az + 2/(1 - z){z 2 - 5 2 ) 

-128f8 + 55 2 U 10 ]ln V , = =}. (39) 

zVA + 5 2 -Az- 2yJ(l - z)(z 2 - 5 2 ) 



" 87r (A 
-{Az, 



z)(z 2 -5 2 ) 



9s 3 5 4 //- 2) 8 (/- 5 2 ) 1 \ (4 + 5 2 -4/ 

x [23045 6 (4 + 5 2 )(144 + 575 2 + 55 4 ) - 11525 6 (6336 + 34245 2 + 5585 4 + 255 6 / 
-1925 4 (12288 - 824965 2 - 391685 4 - 51805 6 - 1755 8 )/ 
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+965 4 (125952 - 1755845 2 - 804085 4 - 88525 s - 2155 8 )/ 

+ 165 2 (73728 - 15797765 2 + 5326405 4 + 3102405 s + 298345 s + 3755 10 )/ 

-85 2 (651264 - 33963525 2 + 2244805 4 + 3239605 s + 248245 s - 2255 10 )/ 

-4(98304 - 24698885 2 + 37414405 4 - 2809285 s - 3272885 s 

-151485 10 + 1055 12 )/ + 2(1179648 - 54927365 2 

+ 11729925 4 - 7960645 s - 2730165 s - 99405 10 + 455 12 )/ 

-(7471104 - 85688325 2 - 22863365 4 - 8642885 s - 1310845 s - 13775 10 )^ 

+ (14909440 - 41123845 2 - 12130565 4 - 332645 s + 11645 8 - 455 10 )/ 

-2(9654272 + 3189765 2 + 1391685 4 + 395245 s + 4475 8 )/ 

+8(1980416 + 2420485 2 + 251205 4 + 8835 s )/ 1 

-64(119296 + 108325 2 + 2455 4 )/ 2 + 1024(1840 + 735 2 )/ 3 - 188416/ 4 ] 

-35 2 (^ - 2) 4 [-1925 s (144 + 575 2 + 55 4 ) + 965 s (1008 + 3045 2 + 155 4 )z 

+ 165 4 (4224 - 63925 2 - 17315 4 - 605 s )/ 

-165 4 (13632 - 9165 2 - 7055 4 - 105 s )/ 

-45 2 (15360 - 564485 2 - 86485 4 + 4335 s + 255 8 )/ 

+25 2 (96256 - 295685 2 - 132805 4 - 3405 s + 155 8 )/ 

+ (16384 - 1935365 2 - 286725 4 + 146805 s + 3995 s )/ 

-5(8192 - 117765 2 - 37125 4 + 6045 s + 35 8 )/ 

+2(10240 + 97285 2 - 27845 4 - 795 s )/ + 512(4 - 475 2 + 25 4 )/ 



„ in zy/A + 5 2 - Az + 2/(1 - z)(z 2 - 5 2 ) 

+256(12 + 195 2 )/ ] In V A =}. 

z^jA + 5 2 -Az- 2y/{l - z){z 2 - 5 2 ) 



9s 3 5 4 z 5 (z-2) 8 (z 2 -5 2 ) 1 \ (A + 5 2 -Az) 

x [23045 s (4 + 5 2 )(48 + 35 2 + 55 4 ) - 11525 6 (2112 + 5765 2 + 1705 4 + 255 s ) 
-1925 4 (6144 - 330885 2 - 72965 4 - 14285 s - 1755 8 )/ 
+965 4 (64512 - 1026885 2 - 162325 4 - 18605 s - 2155 8 )/ 
-165 2 (24576 + 10068485 2 - 6399685 4 - 835845 s - 50065 8 - 3755 10 )/ 
+85 2 (552960 + 30446085 2 - 12879685 4 - 2098645 s - 66085 s + 2255 10 )/ 
+4(98304 - 40919045 2 - 51571205 4 + 28319045 s + 3982165 8 
+28845 10 - 1055 12 )/ - 2(1179648 - 159784965 2 
-51000325 4 + 41536965 s + 3623765 s + 42285 10 - 455 12 )/ 
+ (7471104 - 388300805 2 - 52989445 4 + 32128965 s + 2075885 s - 2855 10 ) 
-(14909440 - 294584325 2 - 31637765 4 + 7699365 s + 158365 s - 455 10 )/ 
+2(9654272 - 59827205 2 - 3887045 4 + 340205 s + 4475 8 )/ 
-8(1980416 - 1936645 2 - 93125 4 + 8835 6 )/ 1 

+64(119296 + 37925 2 + 2455 4 )/ 2 - 1024(1840 + 735 2 )/ 3 + 188416/ 4 ] 
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-35 2 (z - 2) 4 [-1925 6 (48 + 35 2 + 5<5 4 ) + 288<5 6 (112 + 85 2 + 55 A )z 

+ 165 4 (1920 - 37845 2 + 715 4 - Q05 6 )z 2 

-165 4 (6336 - 36605 2 + 289<5 4 - 105 6 )z 3 

-4<5 2 (13312 - 433925 2 + 99285 4 - 8055 6 + 255 8 )^ 4 

+25 2 (92160 - 769285 2 + 69445 4 - 6765 6 + 155 8 )z 5 

+ (16384 - 3225605 2 + 1218565 4 + 109205 6 + 365<5 8 )/ 

-(24576 - 2483205 2 + 968965 4 + 53325 s - 155 8 )/ 

+2(14336 - 409605 2 + 127685 4 + 795 6 )z 8 - 512(36 - 555 2 + 25 4 )z' 
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+ 2 flipped graphs 




Figure 1: Feynman diagrams for e + + e — > 7* — ► Charmonium + cc. 
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Figure 3: Cross sections for e + e~ — > r/ c + cc plotted against the center-of-mass energy. 
Dotted line illustrates the fragmentation calculation and solid line illustrates the complete 
calculation. The cross sections are in units of <j cc = a(e + + e~ — > 7* — > cc) times 10~ 4 . 
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Figure 4: Cross sections for e + e~ — > J/ip + cc plotted against the center-of-mass energy. 
Dotted line illustrates the fragmentation calculation and solid line illustrates the complete 
calculation. 
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Figure 5: Cross sections for e + e" — > Xco + cc plotted against the center-of-mass energy. 
Dotted line illustrates the fragmentation calculation and solid line illustrates the complete 
calculation. 
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Figure 6: Cross sections for e + e" — > Xci + cc plotted against the center-of-mass energy. 
Dotted line illustrates the fragmentation calculation and solid line illustrates the complete 
calculation. 
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Figure 7: Cross sections for e + e" — > Xc2 + cc plotted against the center-of-mass energy. 
Dotted line illustrates the fragmentation calculation and solid line illustrates the complete 
calculation. 
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Figure 8: Differential cross sections of the color-singlet (dotted line) and the sum of color- 
singlet and color-octet (solid line) contributions as functions of the production angle of Xci- 
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Figure 9: Contributions to da(e + e~ — > Xcicc)/dz from color-singlet (dashed line), color- 
octet (dotted line) and the sum of color-singlet and color-octet (solid line) contributions 
plotted against z. 
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